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ABSTRACT 

From radiation magnetohydrodynamic simulations of the solar atmosphere we find a new mechanism for 
the excitation of longitudinal slow modes within magnetic flux concentrations. We find that the convective 
downdrafts in the immediate surroundings of magnetic elements are responsible for the excitation of slow 
modes. The coupling between the external downdraft and the plasma motion internal to the flux concentration 
is mediated by the inertial forces of the downdraft that act on the magnetic flux concentration. These forces, 
in conjunction with the downward movement, pump the internal atmosphere in the downward direction, which 
entails a fast downdraft in the photospheric and chromospheric layers of the magnetic element. Subsequent 
to the transient pumping phase, the atmosphere rebounds, causing a slow mode traveling along the magnetic 
flux concentration in the upward direction. It develops into a shock wave in chromospheric heights, possibly 
capable of producing some kind of dynamic fibril. We propose an observational detection of this process. 
Subject headings: Sun: photosphere — Sun: chromosphere — Sun: oscillations — Sun: surface magnetism — 
magnetohydrodynamics (MHD) 



1. INTRODUCTION 

Filtergrams or spectroheliograms in Ca II H and K reveal 
two main sources of Ca II emission: plages, which spatially 
coincide with the active regions, and the chromospheric net- 
work, whic h outlines the boundarie s of the supergranular ve- 
locity field (Si mon & Leightonlll964 ). Both are also the loca- 
tion of small magnetic flux concentrations in the p hotosphere 
and d ynamic fibrils in the chromosphere (iHansteen et al.l 
120061) . The cl ose relationship between Ca II emission and 
magnetic flux dSkumanich et alj|1975l; ISchriiver et"aniI989l) 
suggests that the magnetic field plays a key role for the 
chromospheric emission. The source of this emission has 
remained elusive, but likely candidates are the dissipation 
of magnetohydrodynamic waves, or the direct dissipation of 
electric currents. With regard to the former process, numer- 
ous studies have been carried o ut based on the app roximation 
of slender flux tubes (see, e.g.. lHasan & Ulmschn eider 2004, 
and references therein). They have greatly expanded our un- 
derstanding of the physics of tube modes, mode coupling, de- 
pendency on the excitation spectrum and the tube geometry, 
shock formation, etc. 

Some of these models consider the generation of tube waves 
by turbulent motions in the convection zone, where an an- 
alytical treatment of turbulence based on the Lighthill-Stein 
theory of s ound generation is used, e.g., in Musielak et al. 
([1981 [1995. 20003 and Musielak & Ulmschneider (2001,) . 
Others consider the generation of transverse waves through 
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the impulse transmitted by granules t o magnetic flux tub es, as, 
e.g., in Choudhuri et al. ( 1993a b), iHasan et al.l (120001) . and 
Cranmer & van Ballegooii en (2005). 

Wave generati on and p ropagation in th ick tub es were 
considered by iZiegler & Uknschneideri (II997al[bl) who 
found the development of a non-axisymmetric longitu- 
dinal body wave and ene rgy leakage by acoustic radia- 
tion. iHuang e t alj (119991) investigated surface and body 
waves of two-di mensi onal magnetic slabs. Mor e recently 
[Rosenthal et al. (2002), Bogdan et al. (2003), Hasan et al 
^005), Hasan & van Ballegooiien (2008), Khomenko et aT 
(2008), and IVigeesh et alj (l2009l) carried out numerical 
simulations of magneto-acoustic wave propagation for an 
environment resembling the network atmosphere including 
small-scale magnetic flux concentrations with internal 
structure. They found that refraction of waves and mode 
conversion may play an important role in the conversion of 
wave energy to heat in the network. 

All the above mentioned models and simulations impose a 
given driving, which is either monochromatic or impulsive, or 
derived from a theoretical spectrum of turbulence, or from ob- 
servations. The focus of the present investigation is the self- 
consistent excitation of a thick magnetic flux concentration 
through the ambient convective motion and specific conse- 
quences of this excitation. 



2. NUMERICAL METHOD AND MODEL CHARACTERISTICS 

The simulations were carried out with the CO^BOLD- 
code dFrevtag et"al] 120021) . This code solves the coupled 
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Fig. 1. — Snapshot of the simulation, showing the full computational do- 
main. Color-scales show the gas temperature. The white solid curves are 
representative magnetic field lines, where the inner pair of field lines indi- 
cate the core region of the magnetic element and the outer pair of field lines 
the inner boundaries of the skin region. The black and white dashed curve 
shows the surface of optical depth unity, t = 1. Arrows indicate the plasma 
velocity. The black dashed box indicates the size of the close-up of Figure|2] 

system of the equations of compressible magnetohydrody- 
namics in an external gravity field and non-local, frequency- 
dependent radiative transfer in one, two, or three spatial 
di mensions. The main featu res of the code are described 
in ISch affenberger et al. (2005, 2006). The two-dimensional 
computational domain extends over a height range of 3 1 60 km 
of which 780km reach above the mean surface of optical 
depth unity. The horizontal extension is 1 1200 km. The grid- 
cell size in the horizontal direction is 28 km, in the vertical 
direction it is 12 km in the upper part, continuously increas- 
ing to 30 km through the convection-zone. The lateral bound- 
ary conditions are periodic in all variables, whereas the lower 
boundary is open in the sense that the fluid can freely flow 
in and out of the computational domain under the condition 
of vanishing total mass flux. The specific entropy of the in- 
flowing mass is fixed to a value previously determined so as 
to yield the correct solar radiative flux at the upper boundary. 
The upper boundary is transmitting, viz., du^r ^/dz ~ 0, while 
d Inpgas/dz = —1/Hp and d Inp/dz — —1/Hp with Hp and 
Hp being the local gas pressure and density scale heights, re- 
spectively. 

The simulation starts with a homogeneous, vertical, unipo- 
lar magnetic field of a flux density of 80 G superposed on 
a previously computed, relaxed model of thermal convection. 
The magnetic field is constrained to have vanishing horizontal 
components at the top and bottom boundary but lines of force 
can freely move in the horizontal direction. At the beginning, 
the magnetic field quickly concentrates in the intergranular 
downdrafts of the convective flow. Subsequently, individual 
flux concentrations merge and after w ICQ minutes, the mag- 
netic field concentrates in a single magnetic 'flux sheet' with 
a strength of 2000 to 2600 G at the optical depth unity within 
the flux concentration. This instant corresponds to the time 
t = in our investigation. Subsequently, the magnetic field 
remains concentrated in this single magnetic element for the 
remaining 86 minutes of simulation time. 

This state, however, is not a stationary one — as a conse- 
quence of the interaction with the surrounding convective mo- 
tion, the flux concentration moves laterally, gets distorted, 
and exhibits internal plasma flow. This interaction excites 
magneto-acoustic waves within the flux sheet, in particular 
longitudinal slow modes. In the following we point to a new 
mechanism of wave generation, which was not considered be- 
fore. 

3. AMBIENT DOWNFLOWS AS A SOURCE OF SLOW MODES IN 
MAGNETIC ELEMENTS 
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Fig. 2. — Time series of close-up views of the magnetic element, demon- 
strating the pumping by external downflows. Colors in the upper panels show 
the vertical velocity, colors in the lower panels the magnetic pressure. The 
white solid curves are representative field lines, where the inner pair of field 
lines indicate the core region and the outer pair the inner boundaries of the 
skin region. The dashed curves correspond to the surfaces of optical depth 
unity, T = 1, the dotted curves to surfaces of equal magnetic and gas pres- 
sure, /3 = 1, above which /3 < 1. Negative velocities indicate downflows. 

Figure [T] shows the complete computational domain with 
the magnetic flux concentration (white magnetic field lines) 
in the middle, the temperature field (color scale), the velocity 
field (arrows), and the continuum optical depth unity (dashed 
contour) for an arbitrary instant in time. The snapshot shows 
warm granular upwellings framed by narrower cool intergran- 
ular downdrafts in the convection zone and a few shock waves 
in the top part of the atmosphere. It also shows two narrow 
downflow channels in the close vicinity outside of the mag- 
netic flux concentration around and below the surface of op- 
tical depth unity. These 'downflow jets' are a consequence of 
a baroclinic flow impinging on the magnetic flux concentra- 
tion from the lateral directions, driven by the r adiative cool- 
ing at the 'hot walls' of t he flux concentration (iDeinzer et al.l 
ll984l;ISteiner et al.lll9 98'). Observational evidences for down- 
flows at the edges of small scale mag n etic flux concen t ration s 
were reported b y Rimmele (2004), 'Ishikaw a et al.l (120071) . 
iLanpangen e t al. (2007), and Rezaei et al. (20071). 

In the present simulation, we observe that these downflows 
are far from stationary. They tend to be present most of the 
time but get transiently enhanced, weakened, or interrupted. 
Sometimes, the lateral inflow carries a preexisting regular in- 
tergranular downflow with it. It then merges with the pre- 
existing downflow channel of the flux concentration, which 
results in a particularly strong downflow. Such advection of 
intergranular lanes to the downdrafts close to magnetic flux 
concentrations are reported to exist from observations of fac- 
ulae by De Pontieu et al. (2006). In the following, we identify 
such transients as a source of magneto-acoustic waves within 
the magnetic flux concentration, in particular of longitudinal 
slow modes. 

3.1. The magnetic pumping process 

Figure 12] shows a close-up of a transient downflow and the 
ensuing excitation of a slow wave over a time period of 234 s. 
The colors represent the vertical velocity in the top row and 
the magnetic pressure in the bottom row. The arrows indicate 
the flow speed, the vertically running white curves selected 
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magnetic field lines of the magnetic flux concentration. 

At the beginning of the time sequence, there exists on the 
left side of the flux concentration a vortical flow, which is di- 
rected towards the magnetic flux concentration near optical 
depth unity, turning into a fast downflow immediately adja- 
cent to the flux concentration. A similar inflow and down- 
flow exists on the right hand side of the flux concentration. 
In the course of time, these flow cells drift in the downward 
direction. The flows exert inertial forces on the magnetic flux 
sheet, which can be seen from the sequence of the magnetic 
pressure, pmag- 

Adjacent to where the flow impinges on the flux sheet, the 
magnetic pressure gets enhanced (e.g., at z w —300 km at 
time 13.8 minutes and at z « —500 km at time 15.1 min- 
utes), while further downstream, where the flow speed is high- 
est, and where the flow starts to detach from the flux-sheet 
boundary, the flux sheet expands and the magnetic pressure 
decreases. This situation leads to an inversion in Pmag as a 
function of height, which moves in the downward direction 
together with the surrounding, gradually downward moving 
vortical flow cells. It presses material within the flux sheet 
in the downward direction. More illustratively expressed, 
the surrounding flow is (unidirectionally) kneading the flux 
sheet in the downward direction. This process entails a down- 
flow within the magnetic flux concentration — while there ex- 
ists a general upflow at time t = 12.5 minutes, it turns into 
a downflow, first in the height range where the kneading is 
strongest (around r = 1), later gradually comprising the up- 
per atmospheric layers up to the top of the photosphere at time 
t = 15.1 minutes. The mechanism by which convective flows 
in the surrounding of a magnetic flux tube may induce a unidi- 
rectional fl ow within the tube was earlier recognized and de- 
scribed by iParkerl (11974b . who called this process 'turbulent 
pumping'.' The principle is used in technical devices such as 
the peristaltic pump. Parker (1974) thought this process to be 
responsible for the intensification of magnetic flux concentra- 
tions. Here we find it to be responsible for the excitation of a 
magneto-acoustic slow mode. 

3.2. The generation of the slow mode 

At first, the kneading process transports all material in the 
downward direction, but this process is not efficient enough 
to prevent that some material can escape this trend and start 
to flow upwards in the center of the flux sheet, which is less 
directly subject to the surrounding flow. This upflow starts 
around r = 1 at time t = 15.1 minutes, while in the photo- 
sphere, plasma still flows at high speed in the downward di- 
rection. The high speed downflow collides with the upflowing 
material of higher density, which leads to an upwardly propa- 
gating wave. Since this takes place in the low /3 regime, and 
since the wave is guided along the magnetic field, it is a slow, 
essentially acoustic wave. Further up it steepens into a shock 
front, visible at time t = 16.4 minutes and z sa 200 km. 

Figure [3] shows two space-time diagrams of the event de- 
picted in Figure |2] The left panel shows the vertical velocity 
(colors) and the gas pressure (contours) horizontally averaged 
over a 50 km thick skin region immediately outside of the 
magnetic flux sheet. The right panel shows the vertical veloc- 
ity (colors) and the magnetic pressure (contours), horizontally 
averaged over a core region of the magnetic flux sheet, which 

' Note that nowadays, this term commonly refers to the transport of mag- 
netic fields by turbulence in stratified media, which is a different process. We 
therefore call this process 'magnetic pumping'. 
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Fig. 3. — Space-time diagram of the close-up of Figure |2] between times 
t = 12 minutes and t = 18 minutes. Left: Skin region of the magnetic 
element. Colors refer to the vertical velocity, white contours to the logarithm 
of the gas pressure. Right: Core region of the magnetic element. Colors re- 
fer to the vertical velocity, white contours to the logarithm of the magnetic 
pressure. Negative velocities indicate downflows. Distances between white 
contours are equal in both panels. The gas pressure of the heavy white con- 
tour in the left panel is equal to the magnetic pressure of the heavy white 
contour in the right panel. Note that /3 is not constant and very different for 
these two contours. The solid black curve indicates optical depth unity, the 
black dashed curve indicates the surfaces of equal gas pressure and magnetic 
pressure, (/9)coro = 1- The strong downdraft in the skin region (red region) 
entails the downdraft in the magnetic element with subsequent shock forma- 
tion starting at about {t,z) = (16 min., —50 km). 

encompasses one third of the total magnetic flux. The height 
z = corresponds to the mean optical depth unity outside of 
the flux sheet. 

A downdraft in the flux-sheet surroundings (left panel) 
starts at a height of z w km and propagates downwards 
while strengthening: it also entails some replenishment from 
photospheric layers. The downdraft front is accompanied 
by an enhancement in gas pressure, which later turns into a 
gas pressure deficit as a consequence of Bernoulli's principle. 
Within the flux sheet (right panel), a downdraft ensues, first 
at time t k, 13.3 minutes closely below the (T)coro — 1 level 
at height z w —400 km. This downdraft cannot be replen- 
ished from the lateral direction (unlike in the skin region) but 
is forced to flow along the magnetic field from higher layers. 
Consequently, a fast downflow ensues from z k, —300 km 
upwards. This downflow soon rebounds and turns into a up- 
flow, developping into a shock front, visible from (t, z) = 
(16 min., -50 km) to {t,z) = (16.6 min., +300 km). At 
z = 300 km, the shock front moves with a speed of ap- 
proximately 15 kms^'^. The Alfven speed at this height is 
w 200kms-i. 

FigureSta) shows the space-time diagram of the vertical ve- 
locity for the full time period of 86 minutes. The colors above 
the height level of continuum optical depth unity, (t)^^^.^ — 1, 
refer to the vertical velocity in the core region of the flux con- 
centration. Below (T)corc = 1, they refer to Vz in the skin 
region, z = corresponds to the optical depth unity outside 
of the flux sheet, {t)^^^ = 1. From this diagram it becomes 
evident that some of the most vigorous downflows in the pho- 
tospheric layers of the flux sheet, e.g., at times t — 5.5, 8.8, 
15.1, 21.4, 62.8, 68.8, and 77.5 minutes, are causally related 
to precursory downflows in the skin region of the flux sheet. 
The downflows in the skin region start at around optical depth 
unity. This is also visible in FigureUb), which shows the ver- 
tical velocity in the skin region to the left (blue) and to the 
right (red) of the flux sheet at a height of around z = 0. The 
solid black curve gives the vertical velocity within the flux 
sheet at the same height. Strong downflows and subsequent 
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Fig. 4. — Top panel: Space-time diagram over the full time span of the 
simulation. The black solid curve corresponds to the (T)corc = 1 surface 
in the core of the flux sheet. Colors refer to the vertical velocity in the core 
region of the magnetic element above (T)corc = 1 and in the skin region of 
the magnetic element below {t)^^^^ = 1. Bottom panel: Vertical velocity 
spatially averaged between —500 km < z < km as a function of time. 
The red curve refers to the skin region on the right hand side of the magnetic 
element, the blue curve to the skin region on the left hand side. The black 
curve refers to the core region of the magnetic element. Negative velocities 
indicate downflows. 

shock formation in the photospheric layer of the magnetic flux 
concentration are often preceded by convective downflows in 
the close surroundings outside of the flux concentration. 

From Figure IHb) one can see that most of the time there 
exists a downdraft in the close surroundings of the magnetic 
flux concentration. However, strong downflow events oc- 
cur sporadically only because they are triggered by the un- 
steady convective motion. Such sporadic events are capable 
of exciting a non-linear, longitudinal wave within the mag- 
netic flux concentration as a consequence of the pumping ef- 
fect. In the wake of this impulsive excitation, the atmosphere 
within the flux sheet oscillates at the acoustic c ut-off period 
dRae & Roberts| fr982'; H asan & Kalkofeij[T999h . This is the 
reason for the predominantly periodic pattern in the upper part 
of Figure |3a), despite of the sporadic, intermittent excitation 
by the ambient downflows. 

4. CONCLUSIONS AND DISCUSSIONS 

We have carried out a two-dimensional MHD simula- 
tion of a network magnetic element embedded in a non- 
stationary environment consisting of a 2.38 Mm thick con- 
vection zone layer, and photospheric and chromospheric lay- 
ers up to 0.78 Mm above the optical depth r = 1. In these 
simulations we found a new mechanism for the generation 
of longitudinal slow modes, which was not considered be- 
fore. Sporadic strong downflows in the close surroundings 
of the magnetic element, 'pump' the plasma within the mag- 
netic flux element in the downward direction by the action 
of inertial forces on the magnetic field — a process that was 
first described by Pa rker! ^1974 ) who referred to it as 'turbu- 
lent pumping'. As soon as the transient ambient downflow 
weakens, the pumping comes to a halt and the downflowing 



plasma from the photospheric and chromospheric layers of the 
magnetic element rebounds, which leads to a slow, upwardly 
propagating magneto-acoustic wave. It develops into a shock 
wave in the chromospheric layers of the magnetic element. 
In the wake of a 'pumping event', the flux-tube atmosphere 
tends to oscillate at the acoustic cutoff-period. We identify 
these pumping events as a crucial mechanism for the excita- 
tion of longitudinal slow modes in flux-tube atmospheres, the 
heating in network and plage areas, and possibly the develop- 
ment of some kind of dynamic fibrils. 

A major shortcoming of the present investigation is the re- 
striction to two spatial dimensions, which might favor the si- 
multaneous presence of downflows on both sides of the flux 
sheet. At this point we do not yet know to what extent the 
pumping mechanism carries over to a three-dimensional envi- 
ronment. Also, we cannot exclude that viscous forces (stem- 
ming from the implemented subgrid turbulence model) may 
add to some extent to the pumping in our simulation. How- 
ever, the viscosity effect is judged minor as there is a rather 
sharp discontinuity in velocity across the (3 — 1 contour, as 
can be seen from Figure |2] 

For a future observational detection of the magnetic pump- 
ing mechanism, we propose to measure Doppler velocities 
within and in the close surroundings of magnetic elements at 
highest possibl e spatial resolution, and fo r different spectral 
lines, s imilar to lLangangen et al.l (120071) or lBellot Rubio et all 
(120011) . Different from these observations however, the mea- 
surements should preferably result in bi-dimensional Doppler 
maps, which would enable to keep track of the magnetic ele- 
ment. They should be supplemented with G-band and Ca II 
H/K filtergrams for the identification of magnetic elements 
and the detection of the chromospheric response, respectively. 
In a time series of such measurements, we would expect a 
transient downdraft in the immediate surroundings of a mag- 
netic element to entail a downflow in the photospheric layers 
within the magnetic element, followed by a fast upflow and a 
brightening in Ca II H/K as a consequence of the shock wave. 
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